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ABSTRACT 



An analysis has been conducted of the capabilities of high 
pressure ratio^ single-stage impulse turbines. Some recent experi- 
mental work has shown considerable promise for these machines. This 
work provides a tentative basis for better understanding the potential- 
ities and peculiarities of such high performance turbines. 

The present analysis indicates that well designed supersonic- 

rotor impulse turbines can be expected to produce static efficiencies 

closely approximating those predicted by well-known (subsonic) curves 

of static efficiency versus ^ • 

isentropic jet speed 

Apparently, there presently exists no generally applicable 
method of blade row loss estimation in the case of supersonic flow. 

An extension to the Ainley-Mathieson loss prediction procedures to 
the supersonic case is proposed. 

Experimental and estimated supersonic -rotor loss coefficients 
of the present report differ little in magnitude from well-established 
results for subsonic blades of similar deflections. 

The effect of rotor incidence, in supersonic flow, has a 
profound influence on pressure and velocity distributions throughout 
the stage. The effect is inadequately explained by a recently proposed 
theory. The present analysis proposes an alternative method for esti- 
mating the effect. 

A logical and straightforward procedure is developed for 
selection of the ’’optimum" design point for a single-stage impulse 
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turbine^ given a required RPM^ pover output^ and rotor blade root 
stress . 

The preliminary design of a large rocket-turbine is con- 
ducted as an exercise in application of the foregoing developments. 

This study was conducted as partial fulfillment of the 
requirements for the Degree of Aeronautical and Astronautical Enginee 
at the University of Michigan. 
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I. INTRODUCTION 



In recent years ^ considerable interest has centered around 
the development of high pressure ratio turbine stages for rocket- 
pump drive and other specialized applications. Desirable character- 
istics are reliability and light weighty with high efficiency as a 
desirable objective of lesser importance. One type of machine attractive 
for these applications is the single-stage impulse turbine. This turbine 
combines low axial thrust on the rotor with potential for high pressure 
ratios and specific work output. The high pressure ratios normally 
result in supersonic flow through the rotor^ and very few experimental 
results concerning these stages are available. 

The present analysis is undertaken to explore the potential- 
ities of such stages with the benefit of some recent test results. The 
study is directed along three parallel^ but distinct topics. These are: 

(1) analysis of flow patterns and loss distribution within the stage 

(2) development of a scheme for selecting the design point to best 
achieve specified design goals ^ and (5) a typical design problem 
intended to illustrate the results of topics (l) and (2). 
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II. ANALYSIS OF LOSSES AND FLOW PATTERNS IN SUPERSONIC TURBINE STAGES 



1 . General Discussion 

While a wide body of theoretical and experimental data are 
available (e.g. Reference l) for predicting flow patterns and blade 
losses at low Mach numbers^ little information is available at super- 
sonic velocities. As a consequence^ turbines designed to operate in 
these regimes frequently fail to achieve the expected performance and 
pressure distributions. Two examples are the turbines of Reference 2 
and 3 wherein the flow patterns^ pressure distributions^ and efficiencies 
differed considerably from estimated design values. In both cases^ the 
results would indicate efficiency was lower than necessary^ had the 
losses been accurately predicted in advance and compensated for in 
the design geometry. 

The method of Reference k was utilized by the N.A.S.A. to 
successfully correlate rotor blade momentum thickness and "specific 
blade loss" with "blade surface diffusion parameter" for one family 
of several transonic-rotor turbines. This method does not include 
the effect of total rotor deflection angle ^ and has not been satis- 
factorily generalized. The method predicts a grossly optimistic rotor 
velocity coefficient for the turbine of Reference 2. The rotor deflec- 
tion angle is l40^ in the latter rotor as compared to only 85 *^ in the 
N.A.S.A. family^ explaining the deviation. 

While the loss prediction methods of Reference 1 are compre- 
hensive^ they are only recommended by the authors in cases wherein 
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the flow is everywhere subsonic. Hence^ an extension in generality 
to include supersonic rotor flow is essential to provide for accurate 
flow channel layout and a more precise estimate of performance. 

Figure 1 presents the performance of the two supersonic 
rotor turbines of Reference 2 and together with that of the N.A.S.A. 
transonic family of References 5 through 10. For comparison^ there 
are the theoretical curves from References 11 and 12 and an unpublished 
curve from the Aerojet -General Corp. The theoretical curves presumably 
are intended to apply to subsonic stages only as no mention is made of 
Mach number levels. Slight deviations between the curves are explained 
by minor differences in the geometry and loss assumptions used in 
deriving them. They may be assumed to fairly approximate the range 
of efficiencies to be expected for single stage impulse turbines 
having negligible Mach number effects. 

Comparison of these curves with the experimental results of 
the high performance turbines shows agreement which is enlightening^ 
if not surprising. It seems that no appreciable deterioration in 
performance may be attributed to Mach effects. The mediocre effic- 
iencies of the turbines of References 2 and 5 are evidently due to very 
high leaving losses^ while the N.A.S.A. family has low leaving losses. 

The turbine of Reference 2 shows performance that is quite 
encouraging when compared to theoretical values. It is no surprise 
that this is the most recent and probably the most carefully designed 
of the three. This turbine is designed for high (l40^) rotor deflec- 
tion and maximum static efficiency. On the other hand^ the N.A.S.A. 
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family of turbojet research turbines are designed for high total 
efficiencies^ which average about 85^* 

The relatively inferior performance of the turbine of 
Reference J appears to result from causes other than the high Mach 
number levels although insufficient data are reported to provide for 
a precise analysis of losses. The performance evidently is degraded 
by extremely severe flow separation from the rotor suction surface. 
This explanation is substantiated by the almost negligible change: 
in flow characteristics which resulted from a major modification of 
the rotor profile. The modification consisted of reducing the blade 
thickness about 40^ by removing metal from the suction surface. The 
negligible resulting flow change indicates the flow was separated in 
both cases ^ and demonstrates the ineffectiveness of profile modifi- 
cations downstream of the separation point. Since the profile 
violates the criteria of Reference IJ^ as discussed in II-5> consider- 
able separation is to be expected. In addition to the separation 
problem^ the design apparently suffers somewhat from relatively 
large tip clearance and "lap". "Lap" refers to a sudden annulus 
enlargement immediately downstream of the nozzle exit. 

Returning to considerations of the turbine of Reference 2^ 
Figure 1 also illustrates the margin of improvement - about two 
percentage points in static efficiency - obtained by increasing 
rotor flare enough to compensate for actual rotor losses. This 
flare increase allowed design impulse conditions to be closely 
approached^ and illustrates the effect of accurate loss estimates 
upon geometry and performance. Comparison of theoretical performance 
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with the best performance of this turbine strongly indicates that 
any losses associated with supersonic effects must indeed be 
negligible. This deduction diametrically conflicts with standard 
works (e.g. Section of Reference l) which indicate a rapid 
increase in losses as the rotor flow begins to exceed local sonic 
velocity. Apparently a phenomenon exists which is similar to the 
"sound barrier" drag increase of aircraft wings. Evidently^ when 
the flow just exceeds local sonic velocity^ strong and near-normal 
shocks exist which cause separation. Once the flow is fully super- 
sonic^ the flow field appears to stabilize to the extent that separa- 
tion losses need be no more severe than in subsonic rotors of similar 
geometry. The design criteria for minimizing rotor flow separation 
are considered in II- 5 . For the present^ it seems reasonable to assume 
that the combined profile^ secondary^ and tip losses of the supersonic 
rotor far outweigh any losses credited to supersonic effects. 

2. Development of a Method for Loss Prediction in Supersonic Turbine 
Blading 

In view of the foregoing examination of Mach number effects^ 
the loss prediction methods of Reference 1 bear re-examination. The 
basic Ainley-Mathieson loss parameter is 

Y = (Total Inlet Pressure) - (Total Exit Pressure) 

(Total Exit Pressure) - (Static Exit Pressure) 

The total pressures are referred to relative and absolute velocities in 
the rotor and stator cases respectively. Following Reference 1^ the 
total loss coefficient for a blade row^ is the sum of profile losses^ 
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Yp^ secondary losses^ ^SEC*^ losses^ Profile losses for 

stator or rotor blading of arbitrary reaction are empirical functions 
of solidity^ deflection angle ^ and thickness ratio. Profile losses 
are easily determined for blades of arbitrary reaction using the 
methods of Reference 1 with Figures 2a and 2b^ which are reproductions 
of Figures 4a and 4b of Reference 1. In the present case attention is 
directed only to the special cases of nozzles and impulse rotor blades. 
In this case the value Yp is determined from Figures 2a or 2b^ as 
appropriate^ for standard thickness — = .20. The profile loss estimate 
for the nozzle or impulse rotor blade in question is then simply 



Yp = 






0.20 



( 1 ) 



The secondary loss equation is developed from relations given 
in Reference 1. From Equation (4) of Reference there is 



(^D^sec /S N 



From Appendix I^ there is Cp = 2(s/c) (tarlP^" tanP^ )cosPjj^ 
and 

(C-r>) = Y (— ) cos^P /cos^P 

^ D^'sec sec^-^ 



m' 



( 2 ) 

(5) 



where the vector mean angle^ 

= tan ^ [| (tan - tan P^)] (5) 



Here p is the gas angle -relative to moving rotor or stationary nozzle 
as appropriate - and x and y refer to inlet and exit conditions 
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respectively. 

Y 

sec 



Combining 2^ 5^ and there results 

^ cos^P^^. ^ ^ \2 

= 4 — I — y (tan P„-tan P ) 
cospia ^ y 



( 6 ) 



The coefficient X is an empirically derived coefficient which accounts 

for the effect of acceleration on secondary losses. Figure I 7 of 

Reference 1 correlates X as a function of hub-to-tip radius ratio 

and (Ay/A^) . Since these data were presumably gather for low sub- 

2 2 

sonic velocities^ it is assumed (Ay/A^) = (W^/Wy) ^ where W is relative 

2 

gas velocity. Based on the foregoing assumption^ (W /w ) should be an 

^ y 

appropriate acceleration parameter independent of M. The assumed 
variation of X with acceleration and hub-to-tip ratio is shown as 
Figure 5. It should be noted that experimental points of the original 
figure presented considerable scatter^ whereas Figure 5 represents a 
best fit of the data. 

The tip clearance loss equation is developed from 



_ 1 Ah /s 

" 2 r/c 

of Reference 1^ page I 5 . Equation (7) is attributed to Carter, and 
was developed on a theoretical basis for small deflections. However, 
Reference 1 points out satisfactory agreement between (7) and tip 
effects at large deflections, hence the extension seems justified 
in the absence of a better theory. Similar to (4) above, there is 
also 



(^d)k 




co^ 

COS^Py 



(8) 
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Upon combining (5)^ (7)^ and ( 8 )^ there results 



, 2(^) p )2 

^ ^ cos ^ 



( 9 ) 



m 



For single stage stators^ P = 0 and (9) simplifies to 



(Y ) = 2(^) 

^ ® h. cos p. 



( 10 ) 



m 



Similarly^ for impulse rotors with p = p = p^ there results 



(Yk)r = 8(|^^-)sin^P (11) 



The secondary loss relation (6) may also he simplified. 
In the case of impulse rotors (6) becomes 



(Y ) = l6x sin^p (12) 

sec R 

In the case of nozzles^ Figure 5 indicates that for any practical 
stator^ Xg = .0055* Accordingly^ (6) becomes 

(Ysec)s = ^(- 0055 )X = .022 (13) 

cosPm cospj^ 

The stator coefficients^ calculated in 11-5^ tend to be 
excessively optimistic if (ij) is used. Furthermore^ the theoretical 
development based on Reference 1 and leading to (ij) has so far made 
no allowance for effects of blade height on the magnitude of stator 
secondary losses. However^ the data of Figure 15b of Reference 1^ as 
well as Reference 9 Q-^d earlier tests ^ confirms that stator secondary 
losses are inversely proportional to blade height^ while solidity^ 
aspect ratio^ and chord have negligible effect on stator losses. In 
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order to resolve the apparent contradition^ (l5) been modified 
by the coefficient 6/h^ chosen empirically based on test results 
presented in II- 5 . Hence ^ the best estimate of stator secondary 
losses is felt to be 



.13 sin^Py 

^ COS 



(14) 



The total losses of a blading row are determined by addition 
of the profile^ secondary^ and tip losses calculated above. That is 



= + Y 



sec 



+ Y, 



K 



( 15 ) 



The Ainley-Mathieson loss parameter may be converted to a pressure 
coefficient or velocity coefficient to simplify subsequent calcula- 
tions . 



One may write 



Y 



p fr p ff p m/p tf -| 

■^x ~ / -^y " 

P ” - P 1 - P /P " 

y y y y 



( 16 ) 



(note P" becomes P' in the case of stationary blading). Then 



= P "/p ' 

y X 



1 + Y[1 






The pressure coefficient is seen to be a function of the relative exit 
Mach number^ as well as the loss coefficient and gas properties. 

The conversion to a velocity coefficient is made using 



^ V (,ti - 



2 /I 







iim 
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which is developed in Appendix III of Reference ik. Solving for 
there results 

t - [1 + (1 - Op^)]' 5 ( 18 ) 

(7-l)My2 

Equations (17) and (18) were programmed on a computer 

and the results tabulated in Table I for 7 = 1.40 and 7 = 1.25« The 
results are also plotted for 7 = 1.40 and appear as Figures 4 and 5- 

5 . Comparison of Loss Prediction Theory With Available Test Results 
The loss prediction theory of II-2 was compared with test 
results for 7 rotors and 5 stators of experimental turbines. These 
represent the currently available results that have been published 
with sufficient information to provide a comparison between predicted 
and experimental loss coefficients for blading rows. The comparisons 
covered a moderately wide spectrum of degree of reaction^ 

deflection angle^ solidity^ aspect ratio and other parameters. While 
the volume of test data is quite limited^ agreement between theory 
and experiment is reasonably good. The average deviation between 
predicted and experimentally determined coefficients is about one 
percentage point. The precision of both predicted and experimental 
coefficients is probably no better than plus or minus one pointy for 
reasons discussed below. The comparative results are shown in Table II 
It was necessary to apply certain assumptions and approxi- 
mations in compiling the data of Table II. The more important of these 
with reasons therefore^ are: 
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(a) The N.A.S.A. rotors of Reference 5 through 9 
were assumed to have tip clearance Ah/h = .02 
or Ah = .042”, since clearance was not specified. 
The errors should he small, as ^ in all 

cases . 

(h) Geometric and velocity characteristics are deter- 
mined for mean radius conditions, deemed repre- 
sentative of the stage. This assumes that three- 
dimensional effects are minor. 

(c ) The effects of Reynolds Number variation are 

neglected. Reference 1 assumes = 2 x 10^, 

while the experimental turbines operated at (Re)^ 
only 2 to 3 times this value. 

(d) The rotor pressure coefficients were not quoted 
explicitly in References 5 through 9* The co- 
efficients were deduced from 




which is developed in Appendix B. While this relation 

is exact, some interpolation was required to estimate 

Hip and from the published turbine maps. The 

stator coefficients were deduced from Reference 9* 

4. Comparison of Typical Rotor Velocity Coefficient Estimates Obtained 
by Different Methods 

Figure 6 displays the estimated variation of rotor velocity 
coefficient, as a function of deflection, for typical impulse blading. 
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The subsonic reference curve is tsiken from Reference 12^ and repre- 
sents a mean obtained from a large body of experimental results. 

The two supersonic curves were obtained following the method of 
11-2^ for a typical geometry. The comparison illustrates several 
interesting features. The magnitude of the rotor coefficient estimate 
differs relatively little between the subsonic and supersonic curves. 
This is not too surprising in view of the comparisons shown in II-l. 
However^ the shape of the supersonic curves at high deflections appears 
somewhat suspicious^ as subsonic experience would anticipate a more 
rapid deterioration of velocity coefficient. Since there is insuffic- 
ient test data to confirm or deny the suspicions^ one would probably 
be well advised to use the more conservative estimate^ or perhaps 
a mean. One also notes the increase of with for fixed losses^ 

Y. This trend is also seen in Figure 5 * 

5 • Detailed Analysis of the Effect of Rotor Profile Geometry on Losses 
Reference IJ reports extensive two-dimensional cascade tests 
of rotor profiles. These tests compared conventional constant- 
curvature steam turbine profiles to several experimental profiles 
designed to minimize flow separation in supersonic rotor channels. 

The inlet Mach number was I .9 and deflection angle was lUO^. 

Schlieren photographs and total pressure measurements indicated 
that the resulting ^^shockless" profiles essentially eliminated 
shock-induced separation from the suction surface. The over-all 
velocity coefficient correspondingly increased from .929 to .952. 
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Reference 2 compared the conventional and shockless profiles 
in a test turbine ^ which also operated with several combinations of 
rotor ”lap" and flare. Exact comparison between cascade and turbine 
tests is difficult due to the limited number of configurations reported. 
Furthermore^ the rotor inlet Mach number was about 1.6 instead of the 
design value of 1.9- The reduced Mach number resulted from off design 
nozzle pressure ratio^ the reasons for which are discussed in Section 
II-6. 

The maximum rotor velocity coefficient achieved in the turbine 
was . 872 . This peak performance was obtained using zero lap^ but 
sufficient flare to approximately compensate for actual losses and 
approach the design exit flow conditions. In another test^ the 
shockless profile increased total head efficiency about 2-l/2 percent- 
age points^ under otherwise identical conditions^ while static efficiency 
changed little if any. These results will be discussed in more detail 
presently. 

The shortcomings of the conventional^ constant -curvature 
passages are thoroughly explored on two-dimensional theoretical grounds 
in Reference I5. Essentially^ it is demonstrated that constant-curva- 
ture passages cannot transform a uniform incident flow into a uniform 
curved flow. Rather^ from considerations of a characteristics net^ it 
is evident that the pressure on either surface undergoes a continuous 
zig-zag pattern of abrupt increases and decreases throughout the length 
of the curved channel. The pressure rises are expected to react with 
the thickening boundary layer and produce separation from the suction 
surface. Cascade tests demonstrated the predicted separation for the 



conventional profiles. 
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A second and complementary fault of the conventional profiles 
is the likelihood of excessive channel width caused by too large a 
spacing between blades. In supersonic rotor channels^ the wider the 
channel^ the larger becomes the uncompensated convex surface acceler- 
ation before the compression waves become incident. Hence^ an 
excessively wide spacing produces an increasingly severe pattern 
of pressure fluctuations. Schlieren photographs also showed a 
definite tendency toward focusing of the compression ”fan" if the 
spacing were too wide. The abrupt pressure rise triggered separation 
and flow breakdown. It should be noted that supersonic rotor profiles 
are determined from channel flow considerations^ and solidity is not 
an appropriate design variable. 

The ^^shockless" profiles of Reference I 5 were designed to 
eliminate the foregoing objections. Their success in cascade tests 
was mentioned above. An example design is discussed and constructed 
in Section IV-5 and Figure 15 . Essential features of the design scheme 
are : 

(a) Dual curvature^ having an entry or transition 
curvature half as great as the central or mid- 
chord curvature. The transition curvature 
extends to a position corresponding to intersections 
of the first characteristics originating from the 
opposite surface. Throughout the central region a 
uniform free vortex flow is theoretically established 
in planes parallel to the machine axis. A similar 



transition at the channel exit reconverts the flow 
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to uniform parallel exit conditions. 

(b) A sufficiently small channel width to reduce to 
acceptable levels the initial^ uncompensated 
acceleration along the convex surface. The 
maximum recommended width is that prescribed 
by a transition curvature having 



(7^), _ 



/ • -1 ^ ^ 
> cos (sin — ) 



( 1 ) 



(c) A channel width sufficiently small to preclude 

deceleration to sonic velocity along the concave 
surface. This requirement is generally less 
stringent than (b)^ unless the incident Mach 
number is near unity. 

A comparison of cascade tests of the two profiles is seen 
in Figure 7^ where the data were collected from Figures 7 19 of 

Reference IJ- Aside from the obvious general superiority of profile 
a closer analysis yields an insight into the nature of the compar- 
ative losses. The general spanwise variation of losses indicates 
the effects of secondary flows ^ wherein the end wall boundary layers 
are "centrifuged" toward mid-span by secondary vortices (see Art. 

15-7 of Reference 12 for complete discussion). A fair estimate of 
purely profile losses is probably obtained from the curve peaks ^ 
near root (and tip). The greatly increased losses in the central 
span region point out the predominance of secondary losses -- in 
qualitative agreement with the estimates of Table II for the rotor 
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of Reference 2. The spanwise uniform superiority of profile B may 
be explained by a sharp decrease in profile losses with secondary 
losses essentially the same in both cases. Hence one concludes 
that the improvement shown by profile B is due to a reduction in 
profile losses caused by virtual elimination of shock-induced 
suction surface separation. 

The ratio of losses between the profiles^ in cascade^ is 

(^AZ^B)cascade " ~ obtained from Figure 5 using the reported 

• 57 

velocity coefficients. From Equation III-2-2, static efficiency is 
directly proportional to (l + for Impulse turbines with all other 
parameters assumed fixed. Thus, the relative improvement in due to 
profile B, if duplicated in the turbine, would be a. factor of only 

„ 1 . 012 . 

1 + .929 

The reported data of Reference 2 allows at least a rough 
comparison of losses obtained using profiles A and B in the test 
turbine. Total efficiencies were .6l0 and .6j5 respectively, with 
only the rotor profiles changed. From Equation II-J-l, (C'p)-g/ (Cp)^ = 
1.06 for stator velocity coefficient and total pressure ratio fixed. 
Also, using Figure i)- for My ~ lA, together with an estimated = .84 
from Table II, one may show ~ Using Figure the 

improvement from, say Y = .84 to Y = .84/1.18 = .71 is equivalent to 
about .013 increase in Because is proportional to (l + 'kp), 

it is understandable that the static efficiency improvement, due to 



the improved rotor profile, was indeed negligible in the test results. 

Several important observations may be made from the results 
of the cascade and turbine tests of the two profiles. It is obvious 
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that the design features of the shockless profile B prove beneficial 
in all cases. The relative improvement of profile B is donsiderably 
less in the turbine^ however^ under this single set of test conditions. 
This result^ together with a velocity coefficient reduction from 
'^Cascade ” *952 to = . 872 ^ may be at least partially explained by: 

(a) the likelihood of moderate flow separation in the 
turbine rotor^ while virtually none occurred in 
cascade tests of profile B 

(b) the degrading effect of tip losses^ not present 
in cascade^ which also are proportionally more 
detrimental to rotor profile B. 

(c) off-design operation of the rotor^ which 
operated at R, =1.6 instead of 1.9- 

(d) other three-dimensional effects^ such as radial 
pressure gradients^ which would be expected to 
degrade the performance of moving blades. 

In addition to the effects of profile changes. Reference 2 
also illustrates another important facet of rotor losses and the 
necessity for their accurate estimation during design. In one test 
series, an improvement from t] = .50 to .52 and = .844 to .872 was 
accomplished by a slight increase in the rotor flare. Both tests 
used the "shockless” profile. In the first case the flare was based 
on an estimated = -91 and, due to the magnitude of the actual 
losses, the exit swirl angle was constrained by continuity to be 
considerably less than its design value. The increased flare approx- 
imated actual losses and produced the increased exit swirl and efficiency. 



-18- 



It is evident that an accurate estimate of the actual velocity 
coefficient may be considerably more important than the question- 
able ability to increase it a few percentage points . 

6- Analysis of the Effects of Rotor Incidence on the Supersonic Flow 
Through Turbine Blading 

On the basis of experiment^ Reference l4 concludes that rotor 
blading with non-zero edge thickness^ in supersonic flow having a sub- 
sonic axial component^ is constrained by continuity to operate at a 
unique incidence relative to the oncoming gas flow. This incidence 
detemines the gas direction relative to the rotor. Likewise^ the 
absolute discharge angle of the stator, as well as effective stator 
exit area, is determined uniquely and is virtually independent of the 
stator geometry. The rotor incidence effect, if uncompensated, can 
drastically reduce the stator discharge area, pressure ratio, and 
Mach number. The stator pressure ratio of the test turbine of 
Reference 2 was reduced to between l/’2 and 2/5 of design value, for 
example, due to existence of rotor incidence of 4.5^- It is evident 
care must be taken to estimate the rotor incidence and thus determine 
the required blading geometry to minimize losses and obtain the 
desired pressure and velocity distributions. 

References 2 and l4 set forth a theory for determining the 
unique incidence, i = P 2 ~ This theory will be stated briefly, 

and its validity examined, in view of the importance of the incidence. 

The incidence is postulated to be that for which the expan- 
sion, corresponding to a Prandtl -Meyer turn of ^2 ~ exactly 

equal to the geometrically available expansion ratio (s cos P^-t)/scos 



There is assumed to be a slender wedge leading edge with attached 



shocks and negligible total pressure loss accompanying the expansion. 
Hence ^ the problem as postulated resolves into simultaneous solution 
of the continuity and Prandtl-Meyer relations. 

The Prandtl -Meyer relation in this case simply v + P = 
constant or 



= v(M„ ) + O 2 -Pj) = v(lvV ) + i 



( 1 ) 



A convenient form of the one -dimensional^ steady continuity equation 
may be developed form 



p VA = ( — ) (M nZ/RT) a = constant 
RT 



Hence ^ 



PMA 

Vt 



= constant 



(2) 



(3) 



But T = T" , and P = 

(1 ♦ M„2)r/r-i 

Since T" is constant through a constant diameter rotor. 



P"MwA 

(1 + j^2)7+1/2(7-1) 



= constant 



(M 



Examining stations (see Figure 9) 2 and 5 , just upstream 
and downstream of the leading edge, one writes 



^ _ V ^ ^7+1/2 (7-1) 



P," ^ 4. Zzi M,,_2 






( 5 ) 
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From the blading geometry^ there is also 



A cos - t/s 

^ - 3 

A 2 cos ^2 

Using Cp = Pj'*/P 2 " and combining (5) and (6)^ there results 



(6) 



t 

S 



cos ^ 7+l/2(7-l) 



Equations (l) and ( 7 ) may be solved for two unknowns. 

Typically^ one knows and t/s and wishes to find and ^ 2 * 

The incidence is thus determined for a specified incident Mach 

number and rotor blade profile. 

Reference l4 illustrates the theory for an example rotor 

having = 70° and = 1-9* Reference 2 compared estimated and 

experimentally deduced incidence^ with fair agreement^ for P^ = 70 ^ 

and = I. 6 I. See Figure 8. Isentropic conditions are assumed 

2 

throughout . 

Upon close scrutiny^ however^ the above theory may be 
shown to become absurd as the rotor blade angle p^^ is reduced 
somewhat. The anomaly is easily demonstrated by a simple numerical 
example . 



Suppose the (relative) incident flow occurs at =1.6. 
The corresponding Prandtl-Meyer angle for 7 = 1.4 is V 2 = l4.9°* 
Further suppose a fixed incidence of any reasonable magnitude say 
5 °. Uow using Equations (l) and ( 7 )^ calculate the allowable edge 
thickness t/s^ for various P^. Assume isentropic conditions^ 



Cp H 1.0. 
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By (l), ) =17-9° which yields a corresponding 



= 1.705- Equation (7) becomes 

t/s = cos P, - cos (P, + 5°) 

_L * 



703 



11^ ^ 

1 + .2(1.703)^ 
1 + .2(1.60)2 



= cos - 1.075 cos Oj + 5°) 



( 8 ) 



For pj = 70 ^^ t/s = .027* However^ upon reducing it is seen that 
t/s decreases monotonically until t/s = 0 for = 52^. The same 
conclusions are reached regardless of the selected numerical values. 

In particular^ it is evident that breakdown occurs at even larger 
if Cp < 1.0 or > 1.6. 

The physical explanation of the breakdown of the theoiy is 
that^ at sufficiently small blade angles^ the turned and accelerated 
Prandtl -Meyer flow (assumed) fails to encounter the requisite geometric 
expansion necessary to satisfy continuity. In addition to the fore- 
going arguments^ it is also evident that the isentropic^ irrotational 
Prandtl-Meyer flow assumptions cannot hold in an actual turbine rotor. 

In article 9*10 of Reference 12^ Vavra has pointed out the inherent 
danger in applying two dimensional analysis to such three dimensional 
flows. 

Aside from the merits of the foregoing theory^ rotor incidence 
remains a physical reality with very important implications. Hence ^ a 
more satisfactory explanation of its nature must be sought out. 

One must agree the continuity relation (7) is valid. Having 
decided the Prandtl-Meyer relation is not applicable in this case^ 



there remains one more unknown than available equations. 
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Expanding the Mach number term of (7 ) in a binomial 
series^ it may be shown that 



Mw. 



1 + ^ = 1 + e - 1) + 9(e2) +. 



J 



Mw:? 

where e = ^ - 1. 

Mw2 

Substituting in (7)^ there results 

+ / A 1 ^2^- 1 ) 

t/s = cos p, - 2 

Cp 



cos Pg 



(9) 



From ( 9 ) it is evident that for moderate incident Mach 

numbers^ the magnitude of the incidence is relatively insensitive to 

small differences between and M^_. 

5 2 

Equation (9) was compared with the experimentally deduced 
incidence of Reference 2. For design speed conditions^ one is able 
to deduce C;^l + 1.9s) = .84. Assuming e > 0^ one then concludes that 



.84 < Cp < 1.00 ; 0 < € < .10 . 



As stated earlier^ a second valid equation in the unknowns 
is apparently lacking. Since it seems probable e « 1 in cases of 
interest^ one might assume e = 0 and rewrite Equation (9) as 

t - cos do 

- cos ± (10) 



The theoretical incidence of Equation (lO) and of Reference 
14 was compared with the experimental incidence of Reference 2. The 
comparison is shown in Figure 8. There is fair agreement^ at the 
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design pointy “between the theoretical estimates for isentropic 
conditions. Furthermore^ Equation (lO) is seen to agree with 
experiment for Cp = .84 at design speed and .90 at zero speed. 

This trend seems reasonable and agrees with the expectations of 
the authors of Reference 2. Higher relative losses at design speed 
are attributed to the reduced Mach number with correspondingly 
stronger^ and possibly detached^ shocks. 

By examining the design speed case further^ one may deduce 
the apparent proportion of total rotor loss_, occurring at the leading 
edge^ which results in the 4.5^ observed incidence. Over-all rotor 
losses were reported as = . 872 . Using Figure 5^ this loss corre- 
s ponds to = 1.12. Similarly^ from Figure 4 the apparent lead- 

ing edge loss^ Cp = .84^ corresponds to (Y-|-p)p = .25* Hence ^ the 
apparent ratio of leading edge to total rotor loss is . 25 / 1*12 = 2 / 9 * 
As noted above ^ this apparent ratio represents an upper limit of 
actual leading edge losses and^ in effect^ incorporates the effect 
of any flow acceleration^ which may occur at the rotor entrance. 
From Equation (lO)^ for a given blade (or gas) angle and incidence^ 
leading edge losses would be expected to increase with t/s. The 
foregoing considerations are utilized to estimate incidence in the 
design problem of IV-5 • 

The result of a single experiment hardly constitutes a 
reliable design rule. However^ there is provided a reasonably con- 
sistent account of the interaction between gas angles^ blade angles^ 
and rotor geometry. As mentioned previously^ care must be taken to 
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accurately estimate incidence and preclude serious rotor-stator 
mismatching with resulting losses. One test cited in Reference 2 
clearly indicates the criticality of the problem. Rotor blades 
having negligible edge thickness were run. During the course of 
the test^ the nozzle pressure ratio decreased from l6/l to Ij/l^ 
corresponding to the accumulation of only a thin carbon deposit 
on the rotor blades. Nevertheless, the buildup caused the 
incidence to evidently increase some 2^ in order to produce the 
observed effect. 



III. A RATIONAL TECHNIQUE FOR DESIGN POINT 
SELECTION FOR SINGLE-STAGE 
IMPULSE TURBINES 



1. General PiscuBsion 

Given a design goal in terms of RPM and power required^ 
expansion ratio, working fluid, and thermodynamic inlet properties, 
one is faced with the task of selecting a design point to satisfy 
the requirements. Literally an infinite number of design points, in 
terms of blade speed and gas angles, can meet the requirements. Hence, 
the designer must choose the most desirable combination of high 
efficiency, compact size and weight, and tolerable stress limits. 

An obvious approach is a trial and error process. This results in 
laborious calculations of debatable final success in determining 
the ’'best" design point. 

A single-stage turbine, especially one with a high expansion 
ratio, will tend to have sizeable leaving losses. Hence, efficiency 
is determined primarily by blade to isentropic jet ratio as seen in 
Figure 1, while the blade speed is limited by material stress con- 
siderations. Experience and the mission requirements generally 
dictate desirable blade materials and allowable stresses. Thus one 
might be well advised to restrict the infinite number of possible 
design points to only those which satisfy some specified blade root 
stress. One could then proceed to decide upon a desirable size, 
weight, and efficiency combination to meet design requirements. 

The following method is one way of logically deciding upon 
the ’’best" design point. It requires neither machine calculations 
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nor infinite human patience. The technique is applicable to single- 
stage impulse turbines^ totally independent of Mach number levels. 

The loss coefficients which enter may be estimated by any means what- 
soever^ and are not related necessarily to loss prediction methods of 
II-2. The nomenclature of the development is defined in Appendix A. 
The method is applied to a design problem in IV, 



2. Static Efficiency as a Function of Velocity Coefficients and Rotor 
Deflection 



In terms of Figure 9^ Euler’s Turbine Equation may be 
expressed as 



T \ = 2 (^2 ^2 



( 1 ) 



For symmetrical impulse blading of the present case^ 



P2 = = P Wg 



Equation (l) may be rewritten. 



T| = 2 JL (1 + tjj) Wj sltl Pj 






or 



T] = 2 — ^ (l + sin 3 



By definition, = Vp/C , and (2) may be expressed 
s ^ o 



( 2 ) 



B = 2il/ (1 + V ) e £ sin 3 



(5) 



From the geometry of Figure 9^ 



0 = n/i - cos 23 - S^sin 3 



(M 
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Hence (5) may be expressed 

•q = 2 ^^ (1 + sin p h{ 4 l - - <Jsin P) (5) 

\|Tj^ is a known function of both P and g ^ in general, tg is inde- 
pendent of g ^ but weakly dependent on P - through OL2* Equation ( 5 ) 
allows calculation of efficiency for selected design point coordinates 
8 ^ p. While Equation ( 5 ) might be considered to "organize" the 
search for a suitable design pointy it places no restriction on the 
blade stresses. 

Incidentally^ Equation ( 5 ) also Illustrates the relative 
importance of tg compared to For example^ a one percentage 

point improvement in to should be roughly equivalent bo a four point 
improvement in since t « f-o ^ 1 * 

■Tv SR 

5. Limiting Stress Requirement 

For a selected blading material of a specified taper ratio^ 
including constant cross section blading as a special case^ one may 
show that blade root stress is proportional to the product of radial 
acceleration and blade height. If bending stresses are neglected^ 

Tt2 

a - ( 6 ) 

R 

By continuity^ the mass flow rate is^ 
o 

W = density x axial velocity x annular area^ or 

W -- cos P Rh (7) 

The exit density^ depends on pressure and temperature at 

the blade exit. The pressure may be assumed independent of the vari- 



ables and equal to exhaust pressure^ provided the axial component of 
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velocity is subsonic. The temperature will vary somewhat with 
“H = Ti If one neglects this variation, Equation (7) becomes 

W ~ cos p Eh (8) 

° 1 

For fixed power output, ^ ~ tJ’* Equation 

(8) becomes 



— ~ ■'1^13 Wo cos p Eh 

Substituting, 

h - i 

T]^j^ W 2 R cos P 

into (6), there results 

U 2 1 

^ ~ Wg cos p 



( 9 ) 



( 10 ) 



One notes that U/R, the angular velocity, is fixed in the 
present case. Further, letting ^ S ^ by assuming 

independent of S and , Equation (lO) becomes 

^ - 1 — 

TitR cos p 



or 



^■n^R S P = B = constant 



( 11 ) 



The parameter B has stress units and may be viewed as an invariant 
property of the particular design specifications and working fluid. 
Experience indicates that B is indeed virtually constant over quite 
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a sizeable range of S ^ P variation^ as illustrated in IV. 

If one now selects a specific working stress b" as a design 
objective^ there results D = — = t]\|/ 5 cos P (l 2 ) 

O' ^ 

Thus^ for a specified stress^ the parameter D is a unique constant. 
Furthermore^ any design point satisfying (12) automatically satisfies 
the imposed stress requirement^ b". 

k. Combination of Performance and Stress Requirements 

By combining the results of Equation (5) and (12)^ 

2 i|/s^(l + tp) sin 3 S ( n/T"^ ^cos^P - Ssin P) = ( 15 ) 

® icos P 

or 

S ( ^yl -g^ cos^p - Ssin P) 2^ sin P cos p (1^4-) 



F' 



F” 



For convenience^ one may define parameters 

^ D 

2^s^ + ^jij) sin P cos P 

= sfl - - Ssin P) 



(l4a) 



(l4b) 



Note that F” is a universal function of S and P and could be tabulated. 

Since S/P) and is assijimed independent of P^ while 

D is a known constant (for a specified w) ^ one may evaluate Equation 
(lU). There remains an unlimited number of admissible S ^ P coordinates 
which are known to satisfy the imposed stress limit. The problem is 
the selection of the most attractive S ^ P that represents a final 



design point decision. 
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As an aid in this decision one might examine some variable 
additional relationships. 

For example^ from (15), 



Tl(S^P) = 



D 



( 15 ) 



Scos p 

The rotor relative Mach number variation may be examined by noting. 



Wp ^ SV2 
*7 /RTg n/ /RTg 

Furthermore, if V2 = ^3*^0 independent of P, as previously assumed, 
it follows that Tg is likewise constant. Hence one notes. 




\ ~ ^ (16) 

2 

Variation in radius is simply determined from 

U ~ R ~ 0 (17) 

Reference 15 gives a useful empirical estimate for single-stage rocket- 
turbine weight as 

m,lb. = TOrt (R,ft.)^ = 1.53 (R,in .)2 (18) 

Finally, since a ~ Rh ~ 0 h, for constant RPM, blade height variation 
may be deduced from 

^ = constant (I9) 

O' 



5. Procedure for Application 

Given power and RPM requirements, together with available 
expansion ratio and turbine inlet conditions for a specified working 



- 51 - 



fluids one may proceed in an orderly manner to select the "best” 
design point, as follows: 

(a) Based on experience and operational objectives, 
one makes a tentative selection of blade material 
and one or more values of allowable centrifugal 
blade root stress. 

(b) a trial design point is selected in terms of an 
arbitrary velocity triangle. As a minimum, one 
evaluates, for the trial design point: 

T) 

O 

W p 

2 

« % 

h 0 




(C) 

(D) 

(E) 



A single calculation of this type should suffice 
provided the calculated stress is within, say + 50 
percent of the selected values, a. 

Calculate the invariant B - Equation (lII-^-ll). 
Calculate D corresponding to selected a - Equation 
(lII-i^-12). 



Plot F’ against P (for a reasonable regime of P) - 
Equation III-4-l4a. Consider constant as calcu- 

o 

lated in step (b). Consider ~ 

using constant as calculated in step (b). The 
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curves F* theNrepresent constant stress lines. 

See Figure 10 for an example. 

(f) Plot F” against for selected S = constant - Equation 
II-4-l4b. Intersections of the constant £ and constant 
a curves determine unlimited compatible ( S ^ a) 
design point possibilities. 

(G) For any desired ( £ ^ a) ^ calculate: 



T) 


Equation (lII-U-15) 


- 


Equation (lII-4-l6) 


9 


Equation (lH-4-4) 


R 


Equation (lH-4-17) 


U 


Equation (lH-4-17) 


m 


Equation (lII-4-l8) 


h 


Equation (111-4-19) 


(H) As an additional 


aid in visualizing the variations^ 



plot T] against P for a = constant. The addition of 
S = constant and 0 = constant points provides for ready 
visualization of the variation in t]^ 
h along the constant stress curve. See Figure 11. 

(l) From the plot of (H)^ one may then select the most 
attractive combination of features to fix the design 
point a). Estimate a-nd h 

for the selected design point. Repeat (b) to verify 
the estimates above. One notes that the approximations 
introduced into the technique are now discarded. A 
second iteration of (C) through (l) would appear unnecessary 
unless (l) indicates intolerable lack of agreement between 



the rough and smooth design point estimates. 



IV. ILLUSTRATIVE DESIGN PROBLEM 



1. Specifications and General Design Philosophy 

An example of a preliminary design is presented to illustrate 
the application of previously discussed methods of loss prediction 
and design point selection. The specifications and gas properties 
arbitrarily chosen are those which might correspond to turbodrive 
requirements of a large liquid propellant rocket using reasonably 
energetic chemical propellants. 

Assumed output requirements are: 28^800 SHP at 15^000 RPM. 

Assumed gas properties are : 

molecular weighty On = 12.0 Ib/mole 

specific heat ratio^ y = 1.25 



constant pressure 
specific heat 



c = 
P 



total inlet pressure P* = 



.827 Btu/lb-°R 
450 psia 



o 

exhaust static pres- 
sure P|^ = 50 psia 

gas constant^ R = 128.8 ft-lb/lb-°R 

total inlet temp- T' = 2200°R 

erature^ ^ 

The turbine is to be of the single-stage^ impulse type to 
minimize axial thrust on the rotor bearings. A reasonable compromise 
between static efficiency^ size and weighty and simplicity is to be 
achieved. 
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2 . Design Point Selection 

The design point selection method described in Section III 
is to be followed. A single^ arbitrary trial design point serves 
as a basis for selection of the final design point. 

The available isentropic energy is^ 

Dt 

' o 

The isentropic jet velocity^ 



^IS - [1 - (p)' *"1 - -827 X 2200 (1 - - 760 



Btu 



C = n/ 2 gJAH^„ = 'j '5 X X 760 = 6160 ft/sec 

O lb . 



The isentropic Mach number, 



^^ 2 ^ IS 



= J - 1 ] = ^/8 [ 15 -'"'' - 1 



] = 2.40 



7-1 Pi, 



The subsequent steps correspond to those of Section III- 5 . 

STEP A 

Suppose consideration be given to blade root stresses of 
50^000^ 57^500^ and 45^000 psia. Most suitable blading materials 
are known to have densities of about O.5O Ih/lnP. 



STEP B 

For the first trial arbitrarily select = 65 ^^ 

From Figure 1 ^ the expected efficiency is .60 < t] < .67. The mean 
blade speed is^ U = *5 ^ 6160 = l 846 ft/sec. The mean blade radius is 

R = ^ X 12 = 

2 jt X 15000 
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(l ) Determination of Stator Velocity Coefficient^ 
The mass flow rate is 



o 

W = 



550 

SHP X 
t)AHi3 



28800 X .707 
•64 X 760 



The annulus area_, 

= 2:!tRh2 = 



and 



o 

W 



o 

W 



Co cosa^ 



h. 



WRT 2 



tg P 2 Cq cos« 2 2jtR 
For an estimated tg « .95* 



^2 = ^o' - 



AH 



-US = 2200 - . 95 ^ = 1570°R 

Cp . 827 



k2 



^2 « 



lb 

sec 



ft -lb 



X 128.8 X 1370 R 



lb 

in 



, 95 X 30 X 6160 



ft 



sec 



cos 65 ^ X 2:rr X l4.1 in 



= 1 . 12 ” 



approximate annulus height. 

The nozzle secondary losses are calculated from Equation (lI-2-l4)^ 



vhere 



Thus 



(^sec )s 



.13 sin^gp 

^2 cos a 
^ m 



= tan“^ (i tan a ) = tan“^ ^ tan 63 ^ ^ i^.yC 



(Y ) 

^ sec'^s 



.13 sln^63° 

1.12 cos 1+7° 



.Ik 
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stator profile losses are estimated from Figure 2a to be (lp)g = -06 . 
Total stator losses^ “ ^^sec^s ^^p^s ** Table I-d, 

for M2 » 2.2, the velocity coefficient estimate is = .97^- 



(2) Determination of some velocity triangle components - See Figure 9 - 
V2 = tgC = .974 X 6160 = 6000 ft/sec 



V. 



ax, 2 



= Vg cos ttg = 6000 cos 65° = 2555 

( Vg Sin Og - U ^ 5^ _ 



ax, 2 



2555 



= 54.7 



w = = ^555 = 1^590 

^ ^2 cos 54.7° 

S = Wg/Vg = 4590/6000 = .752 

0 = u/Vg = 1848/6000 = .508 



^2= ^o ^s 



' - f ^ = 2200 - .974^ = 1528 



.827 



R 



Revised ^ 



6000 



N/7gRTp 72.0 n/1528 



= 2.29 



Mv.0 = 



W. 



72.0 n / t 2 2620 



= 4590 = 1.675 



(5) Determination of rotor velocity coefficient, 

Rotor secondary losses are, using Equation (H-2-12), and 
estimating X = .05 from Figure 5, 



(Ysec)R = sin^p^ = sin^ 54.7° = ^ 



Tip losses^ assuming Z^/h = .05^ from Equation (lI-2-ll) 
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(Y^)^ = 8.0 ^ sin^P^ = 8.0 X .05 sin^ 5 ^. 7 ° = -l6 



K'B. 



Profile losses are estimated from Figure 2 b, ~ 



Estimated total rotor losses are (Ym) 



T^R 



.60. 



With an estimated rotor exit Mach number « 1 . 6 , the estimated 
velocity coefficient from Table I-d is = . 90 ^*-. 

The rotor coefficient compares with an estimated \|/j^ = .898 
from the subsonic reference curve of Figure 6 for P = 5 ^- 7 ° • 

(ij-) Completion of velocity triangle and perfomnance calculations 

Wi^ = ^r^2 = -904 X 1^590 = 5970 
= P 2 = ^^- 7 ° 



''ax .2 - \ ^ ^ 



V,. 



= J 



V, 



+ (W4 cos - U)^ = 2680 



cn, = cos 



, cos-R^) =31.3° 

V. 2680 ^ 



work output \T= = ^2(^ + 1 ^R)sinP 2 ^ l848x 1^590x 1.90i^■ sin 5^.7*^ 

- - 2.5 X 10 ^ 



gJ 



Static efficiency, t) = 



gJ 

V 7 504 

760 



AH- 



= .665 



W = 50k Btu/lb. 



■IS 



This estimate appears reasonable compared with Figure 1 . 



( 5 ) Temperatures^ annulus (rotor blade) heights^ and rotor blade stress 
Since total relative temperature is constant through the rotor 
of constant mear radius^ 
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T " = T ^2^ = Ti, + = Tr" 

2 2gJCp 2gJCp 



Since 



\ ^ V 2 ^ Ti, = + 



(1 - t/) - 1528 . 



2gJCp 



4.15 X 10 



■4 (1 - . 904 ^) 



= 1415 R 
v,,2 



"“*> ■ “ — 



• '4 , 2680 ^ _ O 

* iiJ^ ’ 1T]J7W ’ 

Check work output^ V= Cp AT' = .827 (2200 - 1589 ) = 50 ^ Btu/ll. 

This constitutes a check for consistency of the velocity and temperature 
distributions . 



The mass flow^ 

550 oxrp 

o ^ ^ -7°7 28800 ^ ^o.Ulb/sec 

V 504 



W R To 

The annulus^ or rotor leading edge^ height is hp = Zl — 2 — 

RjtRPpVax^P 

40.4 X 128.8 X 1528 

ho = ^ ^ 7 — T = 1 . 024 " 

^ 2 rt X l 4.1 X 50 X 2555 



li), = ^ ^ ^4 = ho i_ = 1.024 X X —ip = 1 . 207 " 

^"^k^ax.k ^*2+E ^528 .901, 



The blade root stress is calculated assuming a blade of height 
h^ and density . ^0 Ib/in^. The centrifugal stresses only are considered^ 
as bending stresses are expected to represent a small fraction of total 



root stress. 



- 59 - 



The root force^ F = ma = (pA j j3-sec^ ^ ^R_ ^2 f-t/sec^) 

g ft 12 

F = pA, £- 

g 12 6o 

p = F/A = p \ 5_ (50000 = 1910 hl^R psi 

g 12 6o 

for hj^ = 1.207" and R = l4.1" , p = 32600 psi 

(6) Simmary of step (b) results 

Ti = .665 ^ = ,904 

W = 40.4 Ib/sec = 54.7° 

0 = 32600 psi = 1.675 

R = 14.1" = 1.45 

h = 1.207" 0 = .508 

= .97^ S = .752 

STEP C 

B = cTT^tp, S COS 6 -- 8250 psi , the invariant 



STEP D 



D 



( 1 ) 



D 



( 2 ) 




8250 

50000 



.275 



.220 



= .185 



STEP E 

Constant stress lines represented by F' 
are plotted in Figure 10. For this calculation^ 



2^s^'*^r( 1'*' tR)sinPcosP 
.97^ = Constant, and 



tp = 1-^5) is assumed. 



The total rotor losses are 
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esimated following step (B-J) to be 

■ ('^sec + + (^P^R " -7^ sln%. „,t -12 

The velocity coefficient^ is then easily calculated^ as in Step 

(B-5) as a function of variable P with = 1.1^5. A simple working 
plot of constructed from Table I-d is a convenience in 

Steps (e) through (H). 

STEP F 

Curves of for constant are plotted in Figure 10. 

The range .52 < S < .76 represents practical limits of interest and 
increments of .02 are satisfactory. 



STEPS G AW H 



From intersections of constant and o' = constant curves of 

Figure 10^ Equation (III-4-15) was solved for T). Plots of 'r)(P)_ 

cj = const. 

are shown for the selected stresses in Figure 11. While not expected 
to be quantitatively precise^ Figure 11 condenses the -design variables 
into an easily digested presentation. Since ^ h^ depend 

only on 9^ £ ^ the selection of a design point to fulfill a 
specified mission is considerably simplified. 

In the present case_, one can write ^ 



(Eq. II-i^-.l6), 
(Eq. 

(Eq. 111-1^.-17)^ 






9 = plotted for convenience in the 

regime of interest. 



-i+l- 



(Eg. 


111-4-17), 


u « 


1848 


® = 6000 0 






00 

0 


(Eq. 


III- 4 - 18 ), 


m s; 


1.53 


(45.8 0)^ = 3220 9^, lb. 


(Eq. 


111 - 4 - 19 ), 


h. « 


.308 


— X 1.207 = l.l 4 X 10 






9 


^2600 



STEP I 



A working stress must be selected. The blade temperatures 
may be estimated, for a (turbulent) recovery factor of about .9, to 
be. 



T < T + Q ^ 

blade — 2 ’ ^ 2 gJCp 



1528 + .9 X 



4390^ 

4.13 X 10 ^ 



174 T°R = 1287°F. 



From Figure 11 ^ one notes that efficiencies are reasonably good for 
a = 57500^ with relatively little gain to be achieved at higher stresses. 
Experience would indicate a = 57500 to be a logical choice for a 
rocket turbine application of short running time. Ample safety factor 
should be available for a blade material such as Inconel 7O3O at this 
temperature. 

In Figure 11 for a = 57500^ the regime^ .60 < S < .68 appears 
most practical. A rough idea of the trade-off between fuel consumption 
and turbine weight may be quickly made if desired. For example com- 
paring points at = *66 and S2 = shows the latter turbine would 
weigh approximately Zm = 5220(0^ - O^) = 3220 (. 42 ^ - .57^) = 127 lt> . more . 
The corresponding reduction in fuel consumption would be 

^ _ S _ ° . ^^0/7 . 7 8 SHF A _ = .-JO L , x 28800 _ _ 1 _, 

^ ^ tIt Tip 760 .764 .794' 



1.34 Ib/sec . 
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The choice of S= .62 would appear advantageous only, if the burning 
time exceeds 



An 



127 



= 95 sec. 



1.34 

provided the objective is simpipr minimizing the sum of turbine weight 
and total fuel consumed. 

In the light of foregoing considerations a choice was made^ 
selecting: S= p h 60 ^ with the assurance a « 57500. From 

Figure 11 and Step the additional estimates are immediately available 



0 - . 595 

T1 ~ .78 

W ~ 3^ lb/ sec 

1.46 
R ~ 18.1" 
h^ ~ 1.08" 

U » 2370 ft/sec 

m ~ 500 lb. 



The design point selection technique has allowed a logical 
and systematic choice of operating design point. It appears that a 
considerable improvement over the original arbitrary choice has been 
achieved. However^ detailed design calculations must follow to verify 
the choice and remove the inherent approximations . 



Jo Detailed Design Point Calculations 

The absolute gas angle required is 

= cos“^ ( S cos = cos“^(.64 cos 60 ^^) = 71 > 55 ^ 






and 0 = 1.0 sin 0^2 “ S cos ^2 = 71«35^ “ *6^ sin 60° = . 39^ 



(a) Recalculation of stator velocity coefficient^ 

In the manner of Step (B-1)^ must be changed to incorp- 
orate the revised OC2 and nozzle height. For an estimated a .96^ 



To « 2200 - .96^^^ = 1355°R 

.827 



34 X 128.8 X 1333 



30 (.96 X 6160 cos 71.35°) 27t X 18.1 



= .92"^ approximate 



annulus height. For = tan tan 71*35°) = 55.9° 



m 



(Y )s = ^ ^^^ L T1; | ^° = .224 

® .92 cos 55. 9^ 



(Yp) 



P^S 



.06 



Total stator losses, ^^tot^s ” ^2 ^ Table I-d gives 

= >96$> 

It is evident that both (Cp)g and depend on the, as yet 
unknown, value of - estimated to be 2.2. In subsequent calculations 
it is assumed = .965 to determine the actual M2. The value of Cp 
will then be revised, where necessary, to coincide with the newly 
calculated M2. 



(b) Determination of some velocity triangle components 



= .965 X 6160 = 5950 ft/sec \x ,2 " °°® '''^•55° = 1900 

W2 = SV2 = *64 X 5950 = 3810 U = tg Cp© = 5950 X .392 = 2330 



u ^ 2330 
6160 



.378 



a2 = 72.0 -s /1345 = 2640 




3810 

2640 



^2 = 
M2 = 



2200 - .965^ . = 1345°R 

.827 

X2 = 5950 ^ 2.25 
ao 2640 



1.44 



~ 4 i+- 



Estimate 




72.0 nTt4 



.91 X 3810 
72.0 -J 1413 



= 1.28 , with « lhl 3 (Step B- 5 ) 



(c ) Determination of rotor velocity coefficient, '[fp 

(Ygec)p = sin^ 60° = . 36 = 8 ^ -O? sin^ 60° = .18 

(^r)r >12 , and estimated ~ 



From Table I-d, for = 1 . 28 , = .87I . 

An estimated \|r = .881 is obtained for comparison from the subsonic 

li 

reference curve of Figure 6. 



(d) Completion of velocity triangle and performance calculations 



W4 = 



V), = 



,871 X 5810 = 3320 ft/sec 



V , = 

ax , 4 R ax , 2 



.871 X 1900 = 163^ 



n/i 655^ + (5320 cos 60° - 2330)^ a = co 3 ~^ = 18.2° 

1740 

= 1740 



mean radius, R = ^ ^330 x 12 = 17 * 8 " 

2 n X 15000 — 

Work output, . 2330 k 3810 X 1.871 slu 60° _ 

2.5 X 10 ^ 



^ w 575 
^ = TJTrr = = -756 



Mis t 5 o 



o .707 X 28800 

W = = 35-4 Ib/sec 

575 ^ 



The calculated efficiency appears reasonable but perhaps 
slightly optimistic^ compared with Figure 1 . 



(e ) Temperatures and total head efficiency 



* kf ^ fx 0» " = 



2 gJCp 

rp t _ rr , ^4 

■ ^4 + 2iic; 



= i 429 °r 



= 1429 + 



1740 ^ 



4.13 X 10^ 



= 1502 R 
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V- Cp AT' = .827 (2200 - 1502) = 5TT Btu/lb . 

WU 5520 






T2.0nTT), 72.0 n/ 1429 



= 1.22 



M, 



^4 = 



1740 



1740 



72.0 nTT), 72.0 n/ 1429 



= .64 



^4' 



(1 + ^ 50 (1 + .125 X .64^)5'° = 58.5 psia 

W 575 



AH' .827 X 2200 [1 - 
IS 450 



IT = .815 



k. Nozzle Flow Analysis and Profile Layout 

A detailed analysis of the stator losses and required geometric 
layout is necessary. Pertinent assumptions are: 

(a) flow within the nozzles (stations 0-1 ) is considered isentropic^ 

(h) all stator losses occur in and downstream of the exit plane ^ 

between stations 1 and 2^ and a uniform flow is established 

beyond station 2. The losses between 1 and 2 are represented 

by the coefficient ^ = *965? 

s 

(c) continuity is satisfied between 1 and 2, 

(d) tangential momentum is conserved between 1 and 2. 

(e) the flow is adiabatic^ T^^ = \ 

The continuity equation is developed analogous to that 
expressed by Equation (II- 6 - 7 ) except using total temperatures and 
pressures. The resulting relation is 

cos = (|-)g + Cp 

1 -(- .125 



I 

i 

I 



I 



I 



I 
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From previous calculations, Og = 71.55°^ = 2 . 25 . The pressure 

coefficient, 

C s ^ = 30 (1 + .125 X 2 . 25 ^)^*° 



Po' 450 



•773 



Hence, there is 



,t^ ^ (1 + -125 

cos = (i)^ + .0615 ^ 



2\^-5 



Ml 



( 1 ) 



Tangential momentum conservation gives sin o:^ = V2 sin (X^ . 



To' 



\/fgRT^ sin -JT^ 1/ ii sin a^= Mg n/YsR / ^ 



1 + ^ Mn 2 



sina^ 



1 +Z^Mg2 



or 



slnai 1 / 1 Slna, = 






-/ I^ 125 M f 3in 71.35O 



7-1 M 2 n/1 + .125 X 2 . 25 ^ % 



1 + M, 



sin = 1.670 



n/1 + .125 Ml 2 



Mn 



( 2 ) 



For a nozzle with no trailing edge thickness, simultaneous 
solution of (1) and (2) gives Mj^ = 2.150, O.^ = 77*2° . 

For structural integrity a nozzle having 10 ^ "blockage," or 
t/S cos ttg = . 10 , is selected. Equation (l) then becomes 



cos = .0682 



(1 +.125 Mi2)^-5 



Mn 



( 3 ) 



Solution of (2) and (3) gives = 75.3°, M^ = 2.I78 for 
the required nozzle exit conditions consistent with the desired (velocity 
triangle of Figure 12 ) discharge flow. The blade angle turns out to be 
about the largest practically acceptable value. If desired^ the magnitude 
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of could be decreased by providing a slight annular expansion 
between stator and rotor. 

The type of nozzle profile is probably not critical. A 
^'sharp-corner” nozzle of the type discussed in Reference l6 was chosen 
somewhat arbit-rarily . The sharp-corner nozzle^ in general^ provides 
for the shortest possible nozzle length capable of a uniform discharge 
flow. Its properties become much more favorable in cases having larger 
nozzle pressure ratios and lower discharge angles than the present 
example . 

Rather than using a full characteristics network^ the approx- 
imate nozzle contour required was obtained by calculating the required 
throat angle ^ over-all area ratio_, and the approximate orientation of 
the limiting characteristics. The angle of turn at the throat is 
simply i where = Vq^(M 2) the Prandtl-Meyer property angle corre- 
sponding to the exit. 

Vi = tan‘1 ^ - 1) - tan"^ - 1 



\ = 2.178 required angle, 

ratio is obtained from one -dimensional flow theory, 
(e.g. Reference 17 ^ problem 4 . 20 )^ 



The required area 
The relation is^ 



h 

A. 



(^) 7+1/27 



V (^) s(T-) 

’ 7-1 Pi 



TTOT 7 T 



— (-^) 

-1] 7+1 



2 ^7+1/2 (7-1) 



For M-, = 2.178, 



^i'/Pq = (1 + .125 X 2.178 ) 



2 \ 5-0 



10.22 
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and 



(10.22 



,•90 



s/8 (10.22"^'^- 1 



= 2.20 



The schematic nozzle layout is shown in Figure IJ. 

The limiting characteristic of the throat fan^ ab_, stands 



at approximately ^ + p = - 17 - 9 ^ + sin [ 






;] =-17.9° + sin-l(i7^) 



v=l7-9 



19-6 • The limiting reflected characteristic^ bc^ stands at approximately 



p = sin"^( 



2.178 



) = 27 . 4 °. Thus the minimum nozzle length is 



wt 



+ S = (2-8 + 4 . 25 ) = 7.05 wt 



tan 19 . 6 ^ tan 27 . 

However^ due to the required truncation of the nozzle at — - Oiq = l4.7^.» 

the dimension bg must be no less than ^t _ 8.38 w.. Hence, 

tan 14 . 7 ° 

minimum throat to discharge dimension^ parallel to nozzle axis^ is 
eb + bg = 2.8 + 8.4 w.^ = 11 .2 w. ^- . The resulting nozzle profile layout 
is shown in Figure l4a. 

Since the edge thickness^ t^ is l/lO the total widths 

t = 1/9 X 2.20 w.^ = .244 w.^-. Total width is w = 2.444, W-^.. The spacing 

2 . 444 W-J;; 



between blades is s = 



w 



cos aq COS 750 '^ 



9>63 The number of 



blades Z 3 = ^jtR ^ 2 :jt x 17/ Q. = . For a selected throat dimen- 

s 9«o5 

Sion Wq = . 290 "; Zg = 40 blades. The nozzle cascade is illustrated 
in Figure l4b. 

The height of the nozzles (and annulus)^ may be, more pre- 
cisely calculated^ 



o _ 



W R Tp ^ 33.4 X 128.8 X 1343 

o 2 otR 30 X 1900 X 2jt X 17.8 

^ ax y d. 



. 961 ’’ = h. 






The total throat area is ^s^^t = “^0 x . 96 I x .290 = 11. l4 in 
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With the nozzle dimensions fully determined^ the calculated 

stator coefficient of III (a) should he checked. For a thickness ratio 

t/c = .21 and solidity c/s = 1.44^ Figure 2a and Equation II-2-1 gives 
21 

(Yp)s = ^ = .055- The secondary losses are )g 

^ ^ = -272^ where the blade angle rather than the gas 

.961 cos 62.5^ 

angle is employed to provide the most conservative estimate. Total 
losses = .55 exceed the estimated .28 of previous calculations 

Hence ^ from Table I-d^ the estimated = .965 may be about .004 too 

o 

optimistic^ but recalculation of IV-5 seems unnecessary. 



5 • Rotor Flow Analysis and Profile Layout 

A leading edge thickness ratio^ must be selected first 

so that incidence may be determined. For a chord of about 1”^ an 
expected solidity c/s 2.5^ and a desired edge thickness of about 
.020”^ there is t/s = t/c • c/s « .020 x 2.5 = .050. 

The incidence^ may be determined from Equation 

i” C OS So 

(lI-6-lO)^ cos p, = — + — « — ^ . From considerations discussed in 
^ S Op 

11-6^ it is assumed that about one-quarter of the rotor losses occur 
around the leading edge - between stations 2 and 5. 



.66 



165 there results from 



For (Yle)r = - (Yt)r 
T able I-c, (Cp)j^ = Cp(M = l.H, Y = .l65) = ^ 

Hence = cos"^ (.050 + £££_§2_) = 52.8° and the incidence is 7.2° 

^ .90 

The assumed pressure loss has Increased the incidence by about 5.7° 
from the isentropic case. It is of interest that^ were the blade 
angle ^ equal to 60^^ the relative incident gas angle would be 

about 66.1^ and the nozzle pressure ratio would be greatly reduced 



from the desired 15 to 1 value. 
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The rotor profile was laid out following the scheme of 
Reference IJ as summarized in Section II-5 of the present report. 

As discussed in Section 11-6^ the rotor entrance Mach number^ 
assumed equal to = 1.44. It may exceed this value slightly^ as 
suggested by the discussion of II-6. 

By Equation II-5-l^ 






out 






p . I 1 

( — > cos (sln"-^ — — ) = .720 



M ^ ^out - . 2(1 - .720) _ g 



Theoretically a channel width-to-radius ratio of .326 would 
just insure tangency of Mach waves with the convex surface^ while allow- 
ing an excessively large Prandtl -Meyer acceleration on this surface. 



Hence, the authors recommend at least a 20 to 50^ reduction in 

‘ ^^M tran 



/W N 

Accordingly^ a design value ^ ^~^tran ^ selected. For an 



^M 



arbitrary construction scale^ let (^J4)tran ” 1*0^ "^tran 



w. 



cen 



. 22 , 



(ro)tran = 1-H^ ^^i^tran ^ Since the curvature is to be doubled 

in the central section, (rj^)cen = ‘5°^ (^out)cen = (3^i)cen "= -59- 

The profile is constructed, to arbitrary scale of 1.10” 
channel width, in Figure I5. The transition curvature is 15° and the 
leading edge consists of a 10° wedge. This is slightly less than the 
maximum angle consistent with an attached shock. The edges could prob- 
ably be rounded at a slight sacrifice in performance. From construction. 
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the solidity^ c/s = 2.88 and the thickness ratio t/c = ,l 60 . For a 
selected chord of 1 ”^ thickness at the shoulder chamber is approximately 
. 020 ” and the scale of Figure 15 is 5.74/I. 

The number of blades required is calculated from 



2jrR 2 jtR 2jt x 17.8 




522 blades 



The rotor losses may be recalculated more precisely for com- 
parison with = *66 obtained in IV-J. For t/c = ,l 6 and s/c = .35^ 

Figure 2 b and Equation (lI- 2 -l) gives (Yp)p = x .13 = .lOU. The 

^ .20 

revised best estimate of rotor loss is = - 64 . From Table I-d^ 

= tj^(Y = . 6 h, M = 1 . 22 ) = .868 , 



compared to the previous estimate of .87I. From Equation (III-2-5), for 
the final estimates of the final efficiency estimate is t] = .7^8. 
From IV-U, b.2 = = .961", and 



From IV- 2 , a 



H = ^2 7T 
^2 

= 1910 Rh|^ 



1- = ,961 X = 1.172" 

I3U5 .868 

= 1910 X 17.8 X 1.172 = 39800 psi 



B = aT]\lr^ S cos 6 = 39800 X .748 X .868 x .64 cos 6o° = 8270 



Total weight, m r =3220 = 3220 (. 392 )^ = 493 Ih. 



6. Discussion of Resulting Design 

The major numerical results are compiled for comparison with 



estimates of IV- 2 . 
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11 

0 


• 748 


B 


8270 psi 


W 


35-^ Ib/sec. 


P = 


60° 


a = 


59800 psi 


II 

CM 


l.kk 


R 


17.8" 


% = 


1.22 


't's = 


.961 


9 


•392 


't'R = 


.868 


S = 


.64 



The calculated efficiency of . 7^8 agrees well with that 

indicated in Figure 1 for U/C^ = Hence^ a confident prediction 

of efficiencies well above .7 offers encouragement for improvement 

over the previously demonstrated performances of supersonic turbines. 

The Design Point Selection Technique^ leading to Figure 11^ 

proves to be of great assistance in presenting a thumbnail sketch of 

available design choices. The slight variation from the predicted 

values of IV-2 was obviously due to the reductions in both and ^|/ 

n s 

from the approximations assumed. Of particular interest is the observa- 
tion that the "invariant” B changed only .245^ from the value based on 
the original arbitrary design point. 

Perhaps the most questionable area of the problem was that 
of estimating rotor incidence. The theoretical basis is certainly 
open to speculation^ as discussed in II-6. 



V. CONCLUSIONS AND DISCUSSION 



This investigation indicates that well designed supersonic - 
rotor single-stage impulse turbines are capable of rather impressive 
performance. These machines appear to be an attractive compromise 
between simplicity and high efficiency. Utilizing recent test results^ 
an effort is made to better understand the losses and flow peculiarities 
of supersonic stages. Many questions are encountered which cannot be 
decisively resolved without more tests of the quality previously reported 
by the Aeronautical Research Council of Great Britain. 

Analysis of supersonic turbine tests indicate that no appre- 
ciable losses need be attributed solely to Mach number effects^ pro- 
vided certain design criteria are followed. The major rotor blade 
design objective is the prevention of shock-induced flow separation. 
Profiles have been recently developed by the Aeronautical Research 
Council which virtually eliminate separation. In supersonic flow these 
blades have much lower profile losses than the conventional constant- 
curvature blades commonly used in subsonic rotors. Carefully designed 
profiles yield considerably greater improvements in total efficiency 
than in static efficiency. Properly chosen rotor blade flare ^ accurately 
reflecting actual losses^ has a considerable effect on static efficiency. 

Since Mach number level need have no direct effect on blade 
losses^ an extension of the Ainley-Mathieson loss prediction procedures 
to the supersonic case is proposed. Some straightforward modifications 
to the basic procedure are adopted. The extension to supersonic flow 
appears tentatively Justified by the meager amount of available test 
results. 
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The existence of a unique rotor incidence^ in supersonic 
flow through blades of finite edge thickness^ profoundly influences 
pressure and velocity distributions throughout the stage. The 
incidence results from continuity requirements which demand a turning, 
and probably an increase in velocity, at the rotor entrance. The 
explanation of the phenomenon presented in Reference l4 and based on 
Prandtl -Meyer corner flow is shown to be unsatisfactory. An alterative 
explanation is developed based on the assumption of negligible velocity 
increase. More experimental work is required to yield a clear under- 
standing of this important phenomenon. 

A technique is developed to assist the designer in selecting 
the ^optimum*' design point for a single-stage impulse turbine having 
specified RPM, power, and blade root stress. This method is independ- 
ent of Mach number level and the method of estimating blade losses. 

It greatly simplifies the decision by promptly resolving design point 
selection criteria into a readily digested form. 

The preliminary design of a large rocket-turbine is developed 
as a numerical illustration of the foregoing procedures. Blade speeds, 
but not necessarily stresses, are typically very high in order to 
reduce leaving losses. The results indicate static efficiencies 
exceeding .7 may reasonably be anticipated from well designed turbines 



of this type. 
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APPENDIX A 



NOMENCLATURE 



A - area 

B - function defined by Equation III-4-11 

c - chord 

Cp - constant pressure specific heat 

Cp - pressure coefficient; ratio of actual to isentropic total 
pressure 

- isentropic Jet velocity corresponding to stage pressure 
ratio 

D “ function defined by Equation 111-4-12 
F - force 

F^ - function defined by Equation III-4-l4a 

F'* - function defined by Equation III-4-l4b 

g - gravitational constant 

H - enthalpy 

h - blade or annulus height 

Ah/h- rotor tip clearance ratio 

i - rotor incidence angle 

J - mechanical equivalent of heat^ 778 Btu/lb, 

M - Mach number 

m - mass 

N - RPM 

P - pressure 

R - mean blade radius 

R - gas constant 
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blade spacing 
temperature 
blade edge thickness 
mean blade speed 
absolute gas velocity 

gas velocity measured relative to moving rotor 
mass flow rate 
specific work 
channel width 

total inlet pressure - total, exit pressure ^ Ainley-Mathieson 
total exit pressure - static exit pressure 

loss parameter- 
number of blades 

GREEK LETTERS 

Ot - absolute gas angle 

P - relative gas angle 

7 - ratio of specific heats 

s = W2/V2 

g - Mw2 



© s U/V2 

\ - secondary loss acceleration parameter of Figure 5 , 

- wall angle 

^ - velocity coefficient; ratio of actual to isentropic 

velocity 

T\ - efficiency (static) 

a - rotor blade centrifugal root stress 

p - dens ity 



s 

T 

t 

U 

V 

w 

o 

w 

V 

w 

Y 

Z 
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T 


maximum blade thickness 




V 


Prandtl -Meyer angle 




kL - 


sin"^ (^), the Mach angle 




CO 


angular velocity 




SUPERSCRIPTS AND SUBSCRIPTS 




( )' - 


total gas properties referenced to absolute 


gas velocity 


( )" - 


total relative gas properties^ refernc-ed to 
velocity 


relative gas 


( )o - 


stator inlet 




( )l - 


stator exit plane 




( )2 - 


downstream of stator^ after mixing 




( - 


just after rotor leading edge 




( )u - 


downstream of rotor^ after mixing 




( )a - 


profile A of Section II-5 




^ )ax " 


axial component 




( )b - 


profile B of Section II-5 




( )c - 


cross section 




( )cen - 


central or mid-chord rotor channel 




( )i - 


inner or convex rotor profile surface 




( )is - 


isentropic 




( )k - 


rotor tip loss 




^ ^LE " 


rotor leading edge 




( )m - 


mean rotor profile curvature 




( - 


vector mean angle 




^ ^out” 


outer or concave rotor profile surface 





I 

i 
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( ) 
( ) 
( ) 
( ) 
( ) 
( ) 
( ) 
( ) 
( ) 
( ) 



P “ 

R ' 

s ■ 

sec 

T " 

t 

tran 

w ■ 

X 

y 



profile loss 

rotor 

stator 

secondary loss 

total loss^ or total head efficiency 
nozzle throat 

rotor entrance transition channel 
corresponding to relative velocity 
upstream^ or inlet^ position 
downstream^, or exit^ position 



APPENDIX B 



DEVELOPMENT OF EQUATION II-J-l 



By reference to a T-S diagram one may easily show (e.g. 
Appendix B of Ref. U) that 



(Cp)p - 



P 4 ’/Po' 



P2 7Po'(T4'/To'^ 7^ (Cp)g (T4 '/Tq') 



P 4 '/Po' 



,/m 



Ti^'/l'o' expressed in terms of and pij.'/Pq' • 

With Tj^' jg, the temperature corresponding to an isentropic pressure 
decrease top^^', there is Tj^'jg = T^ ' (^4 '/^q ' ^ 



^o’ = ^4; IS 



= T ' [1 - (P4 



,/P 



Since 



T.,.' - T), ' 
Tin = ° ^ 



'T - 



T ' - T, ' 
o 4 



"o’ - "4;is 



To' [1 - 



there is 

Ti^7T^' = (1 - Ti^) + Ti^ (P47 p^') 

and finally^ 



7-1/7 



(cp) 



P^R 



P4VP0' 



(Cp)s [(1-v) + ,iT(Py/Po')(’'-^)/r] r/r-i 
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GAMMA -1.^0 PRESSURE COEFFICIENT AS FUNCTION OF 

Y, M. 7 = 1.40 



-65 




GAMMA =1.^0 VELOCITY COEFFICIENT AS FUNCTION OF 

Y, M. r = l.LO 
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gamma -1.25 PRESSURE COEFFICIENT AS FUNCTION OF 

Y, M. 7 = 1.25 
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GAMMA =1*25 VELOCITY COEFFICIENT AS FUNCTION OF 

Y, M. 7 = 1-25 
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COMPARISON BETWEEN ESTIMATED AND EXPERIMENTALLY DETERMINED LOSS 
COEFFICIENTS OF TRANSONIC AND SUPERSONIC TURBINE BLADING ROWS 
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BLADE VELOCITY ^ JJ_ 

ISENTROPIC JET VELOCITY ’ Co 



Figure 1. Comparison of Theoretical Efficiency Estimates for 

(Subsonic) Single Stage Impulse Turbines with Results 
of Supersonic Turbine Tests, 



PROFILE LOSS COEFFICIENT , Yn PROFILE LOSS COEFFICIENT 



-69- 




(a) 




PITCH / CHORD 



(b) 

Figure 2. Profile Loss Coefficients for Conventional 

Section Blades at Zero Incidence, x/c = .20 
Re = 2 X 1q 5; M < 0.6. (From Ref. l) 
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Figure 5. Secondary Loss Acceleration Parameter. 
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PRESSURE COEFFICIENT, Cp 

Figure 4. Pressure Coefficient as Function of M, Y, 7 = l.U. 
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VELOCITY COEFFICIENT. 

Figure 5* Velocity Coefficient as Function of M, X, 7 = l.if. 
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DEGREES 



Figure 6. Estimate of Rotor Velocity Coefficient - 
Typical Impulse Blade in Supersonic Flow. 
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FRACTION OF SPAN 

Figure 7. Cascade Loss Results of Reference I3. 
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Figure 8. Comparison of Rotor Incidence From Experiment 
and Two Theories - For Turbine of Reference 2 
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stations 








Figure 9* Velocity Triangle Nomenclature. 
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^ -OEG. 



Figure 10 



F' and F” Functions for Turbine Design Problem. 



STATIC EFFICIENCY 
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Figure 11. Constant Stress Efficiency Curves for Turbine 
Design Problem. 
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Flgure 12. Design Velocity Triangle - 
Mean Radius Conditions. 
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’^igure 15. Sharp-Corner Nozzle Geometry Considerations. 

(Not to Scale) 
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Figure l4a. Design Nozzle Profile. 




Figure Xh'b. Shai^-Comer Nozzle Cascade 
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Figure 15. Design Rotor Profile. 
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